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Foreword 


The  Explosive  Excavation  Research  Office  (EERO)  is  embarked  on  a 
program  of  research  in  topical  areas  critical  to  the  overall  technology 
titled  "explosive  excavation."  This  work  was  funded  by  the  Office  of  the 
Chief  of  Engineers  (OCE)  Appropriation  96X3121,  General  Investigations. 
Some  of  these  topical  areas  involve  the  prediction  of  safety-related  effects. 
Effort  is  being  expended  in  these  areas  to  review  all  pertinent  measured 
data  and  all  current  prediction  methods  in  use,  and  to  make  an  attempt  to 
advance  the  state-of-the-art.  This  report  presents  prediction  methods 
that  integrate  and  are  based  on  the  systems  described  in  detail  with  sup¬ 
porting  data  in  two  previous  EERO  reports,  EERO  TR-39,  "A  Revised 
Empirical  Approach  to  Airblast  Prediction,"  and  EERO  TR-40,  "Predic¬ 
tion,  of  Ground-Shock-Induced  Airblast  Overpressures  for  Subsurface 
Explosions  from  Peak  Vertical  Spall  Velocity."  Critical  review  and'com- 
ment  are  invited. 
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Abstract 


Two  methods  are  presented  for  predicting  airblast  overpressures 
from  underground  chemical  explosive  and  nuclear  detonations.  The  first 
method  utilizes  empirical  scaling  techniques  to  predict  both  ground-shock- 
induced  and  gas -vent -induced  overpressures.  Empirical  procedures  are 
also  provided  for  predicting  airblast  from  underground  row-charge  and 
array  detonations  and  from  surface  or  near- surface  bursts.  The  second 
method  uses  peak  vertical  spall  velocities  combined  with  basic  physical 
principles  to  predict  ground-shock- induced  overpressures.  With  the 
use  of  one  or  both  of  these  methods,  it  is  possible  to  predict  airblast 
overpressures  caused  by  detonations  at  scaled  depths  of  burst  from 
zero  to  700  ft/kt1^  in  a  wide  variety  of  media. 
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Nomenclature 


c  =  Sonic  velocity  in  ambient  air  (in  ft/ sec). 

Cq  =  Sonic  velocity  in  air  at  standard  sea-level  conditions  (in  ft/ sec). 

DOB  =  Depth  of  burst  =  distance  from  shot  point  (center  of  charge)  to  sur¬ 

face  ground  zero  -  SGZ  (in  ft). 

dob  =  Scaled  depth  of  burst  to  the  center  of  a  charge  (in  ft/kt1^). 

fmax  =  Maximum  transmission  factor  (may  refer  to  either  ground-shock- 

induced  or  gas-vent-induced  airblast);  fmax  is  used  in  predicting 
airblast  by  the  empirical  method. 

h  =  Height  of  an  elevated  location  above  the  ground  surface  level  (in  ft). 

M  =  Mach  number  of  the  induced  air  shock  at  SGZ  (dimensionless)  (cal¬ 

culated  from  the  peak  vertical  surface  velocity,  Vq). 

m  =  Velocity  profile  constant  (for  all  single-charge  detonations  in  a  given. 

medium)  (dimensionless).  This  constant  is  used  to  relate  the  ground 
surface  spall  velocity  as  a  function  of  inverse  dimensionless  slant 
range  (DOB/S)  to  the  spall  velocity  at  SGZ. 

n  =  Total  number  of  charges  in  a  row  or  array  of  buried  charges. 

Pambient  1  Ambient  atmospheric  pressure  at  the  location  of  an  experiment  -  as 

1  obtained  from  meteorological  or  altitude  data  (in  mbar). 

APg  =  Standard  line  overpressure  as  a  function  of  standard  line  range  (in 

mbar).  Standard  line  points  (Table  1)  are  used  as  reference  points 
in  predicting  airblast  by  the  empirical  method.  For  surface  bursts, 
APS  is  also  used  to  represent  a  standard  surface  burst  overpressure. 

AP  =  Peak  airblast  overpressure  at  any  location  beyond  the  immediate 

vicinity  of  the  detonation;  may  be  used  to  refer  to  either  a  ground- 
shock- induced  or  a  gas -vent- induced  overpressure  (AP  may  be 
expressed  in  psi  or  mbar). 

AP„  1  Ground-shock-induced  peak  local  airblast  overpressure  directly 

above  SGZ  (in  psi). 

AP  (elevated)  =  Peak  airblast  overpressure  at  an  elevated  location  —  above  the  ground 
surface  level  (may  be  expressed  in  psi  or  mbar). 

Rs  1  Standard  line  range  (in  ft/kt1^).  Standard  line  points  (Table  1)  are 

used  as  reference  points  in  predicting  airblast  by  the  empirical 
method.  For  surface  bursts,  Rs  is  also  used  to  represent  a  standard 
surface  burst  range. 

R  =  Range  from  SGZ  to  any  location  at  or  above  the  ground  surface  (in  ft). 

S  =  Slant  range  from  the  shot  point  to  any  location  on  the  ground  surface 

(in  ft). 

SGZ  =  Surface  ground  zero;  the  location  at  the  ground  surface  lying  verti¬ 

cally  over  the  shot  point. 


-vii- 


V 


w 

a 


p0 

pQC0 


e 

i 

II 


=  Peak  vertical  surface  spall  velocity  '  .  .  round  surface  at  any 
location  other  than  SG2  (in  ft /sec), 

=  Peak  vertical  surface  spall  velocity  of  the  ground  surface  at  SGZ 
(in  ft/sec). 

=  Yield  of  a  buried  detonation;  usually  expressed  in  kilotons  (kt). 

=  Dimensionless  airblast  constant  (for  all  single- charge  detonations 
in  a  given  medium);  used  in  the  theoretical  prediction  of  ground- 
shock-induced  airblast;  dimensionless  quantity. 

=  Density  of  air  at  standard  sea-level  conditions. 

-  Acoustic  impedance  of  air  at  standard  sea-level  conditions  (in  psi- 
sec/ft). 

=  Angle  between  the  vertical  direction  and  a  vector  from  SGZ  to  an 
elevated  location  (cos  6  =  h/R);  (in  degrees). 

=  Ground-level  location  in  a  direction  perpendicular  to  the  axis  of  a 
row  of  buried  charges. 

=  Ground-level  location  in  a  direction  off  the  end  of  a  row  of  buried 
charges  (i.  e. ,  along  the  axial  direction  of  the  row). 
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PREDICTION  OF  AIRRLAST  OVERPRESSURES 
FROM  UNDERGROUND  EXPLOSIONS 


Introduction 


PURPOSE 


The  use  of  underground  explosions  in 
construction  applications  is  often  limited 

by  corisiderations  of  damage  to  nearby 

\  ' 

man-made  structure’^  or  to  the  natural  en-  : 
vironment.  There  are  three  major  sources 
of  damage'  from  such  underground  explo¬ 
sions-  material  throw  out,  ground  motion, 
and  airblast.  The  effects  of  airblast  can 
be  a  limiting  factor  at  ranges  greater  than 
several  hundred  feet  from  the  detonation 
site. 

Airblast  from  buried  detonations  has 
been  under  systematic  investigation  for 
two  decades.  Sufficient  data  are  now 
available  to  permit  the  development  of  a 
general  prediction  scheme.  Damage  due 
to  airblast  is  very  closely  correldted 
with  the  peak  overpressure  of  the  airblast 
pulse  (pressure  increase  above  the  local 
ambient  pressure)  for  almost  all  cases 
of  interest.  Therefore,  this  is  the  quan- 

I 

tity  which  .must  be  predicted.  This  re¬ 
port  incorporates  prediction  procedures, 
which  were  developed  in  Refs.  1  and  2  by 
considering  the  data  from  all  previous 
large-yield  chemical  and  nuclear  sub¬ 
surface  detonations,  into  a  general  method 
which  predicts  peak  airblast  overpressure 
as  a  function  of  range  from  the  surface 


ground  zei^o  (SGZ)  of  the  detonation  site. 
The  prediction  method,  is  summarized  and 
complete  instructions  f6r  performing  pre¬ 
dictions  are  given,  together  with  sample 
problems. 

BACKGROUND 

Most  buried  detonations  at  useful  cra¬ 
tering  depths  give  rise  to  t^-o  separate 
airblast  pulses.  The  first  pulse  is  pro¬ 
duced  when  the  shock  wave  from  the  det¬ 
onation  reaches  the  ground-air  interface, 
is  reflected  as  a  tension  wave,  and  forces 
the  overlying  ground  surface  to  move 
upwards.  The  rising  ground  surface 
compresses  the  air  above  it,  producing1  a 
distinct  airblast  pulse.  This  pulse  is  re¬ 
ferred  to  as  the  ground-shock- induced 

\ 

airblast.  The  peak  overpressures  for  the 
ground- shock- induced  pulse  can  be  pre¬ 
dicted  by  using  empirical  techniques  or 
by  using  a  theoretical  method  based  on 
the  theory  of  acoustics  (sound)  and  the  ’ 
predicted  peak  vertical  spall  velocity  of 
the  ground  surface. 

.'Jter  th^  ground  surface  above  the  det¬ 
onation  is  she  ked,  a  mound  of  material 
rises  into  the  air.  'The  mound'grows  for 
a  time,  then  begins  to  disintegrate.  Wherj 
disintegration  occurs,  the  explosion 


\ 
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cavity  gases  beneath  'he  mound  escape\or 
vent  to  the  surrounding  atmosphere.  If 
the  pressure  of  these  gases  at  vent  time 
is  appreciably  greater  than  the  ambient 
pressure  of  the  surrounding  atmosphere, 
a  second  pressure  pulse  will  be  produced. 
This  pulse  is  known  as  the  gas-vent - 
induced  airblast.  At  the  present  time,  its 
peak  overpressures  can  be  predicted  only 
by  using  empirical  techniques. 

In  addition  to  the  two  primary  airblast 
components,  additional  pulses  of  lesser 
amplitude  are  sometimes  observed.  These 
additional  pulses  are  always  smaller  than 
the  dominant  (ground-shock- induced  or 
gas-vent-induced)  pulse,  and  they  may  be 
ignored  for  purposes  of  safety  predictions. 
Figure  1  shows  a  typical  measured 
ressure-vs-time  history  recorded  at  a 
few  hundred  feet  from  a  chemical  explo¬ 
sive  detonation  occurring  near  optimum 
cratering  depth  in  aUuvial  soil.  Note 
that  the  ground-shock- induced  pulse  and 
the  gas-vent-induced  pulse  are  both 
clearly  visible.  The  gas -vent -induced 
pulse  happens  to  be  "dominant"  in  this 
case  (producing  a  larger  peak  overpres¬ 
sure  than  the  ground-^shock-induced  pulse). 

PKKDICTION  MKTHODS 

Most  detonations  buried  at  cratering 
depths  in  unsaturated  media  give  rise'  to 
both  ground-shock- induced  and  gas-vent- 
induced  airblast  pulses.  The  larger  of 
tkese  two  pulses  will  determine  the  max¬ 
imum  peak  overpressures  and  thus  the 
airblast  damage  limitations  for  the  deto¬ 
nation.  Gas-vent-induced  airblast  is  usu¬ 
ally  dominant  for  detonations  at  shallow 
scaled  depths.  However,  the  gas-- vent- 
induced  pulse  i?  very  rapidly  suppressed 
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Fig.  1.  Typical  intermeaiate-range  air- 
blast  overpressure  tracing  (for 
a  TNT  detonation  in  unsaturated 
alluvium),  showing  ground- 
shock-induced  pulse  and  gas- 
vent-induced  pulse. 

with  increasing  depth  of  burst  and  also 

lags  further  behind  the  ground- shock- 

i  \ 

induced  pulse  in  time.  At  great  depths  of  J 

burst,  the  ground- shock- induced  pulse 
becomes  dominant,  and  the  gas- vent- 
induced  pulse  disappears  altogether. 

Since  either  of  the  two  pulses  may  be 
dominant  in  a  particular  case,  it  is  usually 
necessary  to  predict  both  ground-shock- 
induced  and  gas -vent -induced  airblast. 

This  report  presents  two  prediction  meth¬ 
ods.  The  first  method  is  empirical  in 
nature  and  is  based  on  earlier  studies  by 
L.  J.  Vortman.'*  ®  The  empirical  \nethod 
predicts  both  ground-shock-induced  and 
gas-vent-induced  air-blast  for  most  deto¬ 
nations  .n  unsaturated  media.  The  second 

method  is  based  on  a  theoretical  treat- 
6  7 

ment,  and  it  predicts  ground-shock 
induced  airblast  for  certain  types  of  deto¬ 
nations  in  certain  media.  The  theoretical 
method  is  often  more  accurate  than  other 
techniques,  but  it  is  applicable  only  to  a 
limited  number, of  cases. 
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Both  methods  require  the  same  basic 
input  parameters.  First,  the  yield,  W, 
and  the  scaled  depth  of  burst,  dob,  of  the 
charge  must  be  known: 

dob  l2SPth  of  burst  (in  ft) 

|  W  (in  kt)]1//3 

W  =  Yield  of  charge  m  kilotons. 

In  addition,  it  is  necessary  to  specify  the 
type  of  explosive  used,  the  emplacement 
medium,  and  the  approximate  moisture 
content  of  the  medium.  The  theoretical 
method  also  requires  the  peak  vertical 
spall  velocity,  Vq,  of  the  surface  at  SGZ. 

The  empirical  and  theoretical  methoas 
predict  peak  overpressures  as  a  function 
of  range  along  the  ground  surface.  All 
ranges  are  measured  from  SGZ,  the  sur¬ 
face  point  directly  above  the  detonation. 
Empirical  predictions  are  valid  at  or  near 
ground  surface  level  (elevations  less  than 
several  yards).  The  theoretical  method 
likewise  predicts  overpressures  at  or 
near  ground  surface  level,  but  it  can  also 
be  used  to  predict  elevated  overpressures 
(above  the  ground  surface)  if  desired. 

Airblast  can  also  be  predicted  for  rows 
or  arrays  of  buried  charges  which  are  em¬ 
placed  close  together  so  as  to  excavate  a 
223 

continuous  crater.  *  *  The  airblast 
overpressures  for  an  "average"  single 
charge  in  the  row  or  array  are  first  pre¬ 
dicted  by  means  of  either  the  empirical 
method  or  the  theoretical  method.  These 
predicted  overpressures  are  then  multi¬ 
plied  by  a  reinforcement  factor  which 
depends  on  the  number  of  charges  in  the 
row  or  array.  This  procedure  gives  the 
predicted  overpressures  for  the  entire 
row  or  array.  The  reinforcement  factors 
are  determined  by  means  of  an  empirical 
method,  but  may  be  applied  to  either  em¬ 


pirical  or  theoretical  single- charge  pre¬ 
dictions. 

The  correct  choice  of  which  method  to 
use  in  specific  cases  is  discussed  in  the 
following  pages.  The  general  rules  are 
summarized  as  follows:  The  empirical 
method  may  be  used  to  predict  ground- 
shock- induced  and  gas-vent-induced  air- 
alast  from  nuclear  and  chemical  explosive 

detonations  at  cratering  depths  (dob 

1  'X 

**  IS  to  250  or  300  ft/kt  ‘  )  in  unsaturated 
media.  Gas- vent-induced  airblast  is  al¬ 
ways  strongly  dominant  for  detonations 
shallower  than  dob  =  60  ft/kt  ,  and 
ground- shock- induced  airblast  need  not 
be  predicted  in  these  cases.  The  empir¬ 
ical  method  is  also  used  to  predict  rein¬ 
forcement  factors  for  row -charge  and 
array  detonations.  The  theoretical  meth¬ 
od  may  be  used  to  predict  ground-shock- 
induced  airblast  from  nuclear  and  chemical 
explosive  cratering  detonations  in  most 
(but  not  all)  unsaturated  media.  Theoret¬ 
ical  predictions  of  ground -shock -induced 
airblast  are  usually  somewhat  more 
accurate  than  empirical  predictions  when 
the  dob  is  gt eater  than  60  ft '  kt  '  ‘ .  If  the 
theoretical  method  is  used  for  ground  - 
shock-induceo  airblast  predictions  in  these 
cases,  the  gas-vent-indueed  airblast  must 
still  be  predicted  by  the  empirical  method. 
There  are  also  certain  types  of  events  for 
which  ground-shock- induced  airblast  is 
strongly  dominant  and  gas-vent-indueed 
airblast  is  negligible.  The  theoretical 
method  is  generally  used  to  predict  air- 
blast  from  these  events  (although  the  em¬ 
pirical  method  can  be  used  in  some  cases). 
Ground-shock-induced  airblast  is  strongly 
dominant  for  all  mounding  or  contained 

detonations  (700  ft/kt1  ^  >  dob  ^  300  ft/ 
1/3 

kt  '  ),  for  most  well-stemmed  nuclear 
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detonations  of  yield  less  than  5  kt  deeper 
than  dob  =  125  £t/ kt 1  /  ^  in  dry  high-strength 
rock,  and  for  all  detonations  deeper  than 
dob  =  170  ft/kt1/3  in  saturated  media  or 
under  water. 

SURFACE-  BURST  AND  FREE- AIR -BU  RST 
AIRBLAST  PREDICTION 

Airblast  from  detonations  shallower 
than  dob  =  15  ft/kt1/3  >s  very  nearly  equal 
to  the  airblaat  from  a  surface  burst  of  the 
same  yield.  An  approximate  method  for 
predicting  airblast  overpressures  from 
surface  and  near-surface  events  is  in¬ 
cluded  with  the  empirical  prediction 
method.  A  method  for  predicting  airblast 
from  free-air  bursts  is  also  given. 

RANGE  OF  VALIDITY 

Both  of  the  subsurface-burst  prediction 
methods  discussed  here  utilize  straight 
line  log-log  fits  for  predicting  airblast  as 
a  function  of  range.  The  straight  line 
approximation  is  valid  at  intermediate-to- 
long  ranges.  Ground- shock -induced  air- 

blast  predictions  are  accurate  for  scaled 

*  1/3 

ranges  greater  than  600  ft/kt  '  ,  out  to 

true  ranges  of  at  least  several  miles  or 
more.  Gas-vent-induced  airblast  pre¬ 
dictions  are  accurate  for  scaled  ranges 
greater  than  3000  ft/kt1/'*,  out  to  true 
ranges  of  at  least  several  miles  or  more. 

Scaled  range 

(ft/kt1/3)  =  actual  range  (in  ft) . 

IW(inkt)]1/3 


The  true  ranges  at  which  the  predictions 
are  accurate  for  a  typical  chemical  ex¬ 
plosive  detonation  (yield  =  1.0  ton)  would 
be  60  ft  to  several  miles  (ground- shock- 
induced  airblast)  and  300  ft  to  several 
miles  (gas -vent- induced  airblast).  Local 
meteorological  effects  may  become  impor¬ 
tant  at  ranges  beyond  several  miles.* 
These  long  ranges  are  usually  of  no  inter¬ 
est  for  buried  chemical  explosive  detona¬ 
tions  smaller  than  1.0  kt  because  the 
airblast  overpressures  are  not  great 
enough  to  cause  damage. 

Within  the  range  intervals  specified 
above,  predicted  airblast  overpressures 
are  normally  accurate  to  about  ±  40%. 
The  accuracy  may  be  somewhat  lower  for 
media  or  explosives  different  from  those 
specified  in  the  prediction  methods,  for 
highly  inhomogeneous  media,  for  indivi¬ 
dual  charges  of  very  small  yield  ( less 
than  one  ton),  and  for  row -charge  or  ar¬ 
ray  detonations.  Predictions  are  not  ex¬ 
tremely  accurate  for  very  long  or  complex 
rows  of  many  charges. 

The  surface-burst  and  free-air-burst 
overpressure  predictions  are  valid  at 
scaied  „  unges  greater  than  30  to  100 
ft/kt  1  ,  out  to  true  ranges  of  several 
miles  or  more .  The  predicted  over¬ 
pressures  are  usually  correct  within  ±40% 
for  the  explosives  specified  in  the  predic¬ 
tion  method. 


Atmospheric  refraction  and  focusing 
effects,  which  become  important  at  very 
1c  ;g  ranges  for  large-yield  detonations, 
are  discussed  in  Ref.  9. 
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Empirical  Prediction  Method 


PREDICTION  PROCEDURE  FOR  SINGLE¬ 
CHARGE  SUBSURFACE  BURSTS 

Airblast  overpressur  es  from  buried 

single-charge  detonations  at  cratering 

depths  (15  <  dob  <  250  or  300  ft/kW3  ) 

in  unsaturated  media  may  be  predicted  by 

the  method  described  in  this  section.  The 

method  involves  a  scaling  process  whereby 

"standard  line"  overpressure-vs-range 

points  are  scaled  to  the  yield  and  ambient 

atmospheric  pressure  of  the  experiment 

to  be  predicted.  The  scaled  overpressures 

are  then  multiplied  by  a  ground- shock 

transmission  factor,  f  ,  (ground- shock), 

max 

to  predict  the  ground-shock-induced  air- 
blast,  and  by  a  gas-vent  transmission 
factor,  f  (gas -vent),  to  predict  the 

gas- vent- induced  airblast.  These  trans¬ 
mission  factors  are  empirically  deter  - 
mined  corrections  which  convert  the 
scaled  overpressures  to  true  predicted 
overpressures  for  the  detonation  in  ques¬ 
tion.  The  following  are  the  steps  of  the 
empirical  prediction  procedure. 


Step  1 

Calculate  the  scaled  depth  of  burst, 
dob,  of  the  detonation  to  be  predicted: 


dob 


DOB 


(1) 


DOB  =  Depth  of  burst  to  center  of 
charge  (in  ft) 

W  =  Yield  of  charge  (in  kt). 

Scaled  depth  of  burst,  dob,  is  expressed 
in  units  of  ft/kt^3. 


dieted.  First,  select  several  standard 
line  overpressure-vs-range  points  from 
Table  1.  Then  scale  the  standard  line 
overpressures  and  ranges  according  to 
the  following  equations: 


AP 

R 


=  R 


ambient) 

1000  mbarJ 

W  1000  mbar1 


1/3 


1.0  kt  P 


ambient 


(2) 

(3) 


where 

AP  =  overpressure  (in  mbar) 

R  =  range  (in  ft) 

AP  =  standard  line  overpressure 
(in  mbar)  from  Table  1 


K  -  standard  line  range  ( in  ft/ 
ktl/3)  from  Table  1 

^ambient  r  ambient  atmospheric  pres¬ 
sure  at  location  of  experiment 
(in  mbar) 

W  yield  of  charge  (in  kt). 


The  ambient  atmospheric  pressure  , 

^ambient  at  Pr°ject  location  may  be 
estimated  either  from  meteorological  data 
or  from  the  altitude  of  the  location  above 
sea  level.  Normally,  Pambient  =  1013 
mbar  at  sea  level,  and  pambjent  * 
mbar  at  5000  ft  above  sea  level  (a  34 -mbar 
decrease  for  each  1000  ft  above  sea  level). 

In  order  to  insure  the  accuracy  of  the 
predictions,  it  is  necessary  to  scale  three 
or  more  standard  line  points  (R^,  APC). 
These  points  should  be  selected  at  a  vari¬ 
ety  of  Rg  ranges.  The  scaling  process 
will  then  provide  a  set  of  points  (R,  AP) 
which  are  used  in  making  predictions. 


Step  2 

Calculate  several  overpressure-vs- 
range  points  for  the  detonation  to  be  pre- 


Step  3 

Determine  the  transmission  factors 
f  (ground-shock)  and  f  (gas-vent). 

mav  mav  ° 
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Table  1.  Approximate  values  of  standard 
line  overpressure  APa  as  a 
function  of  range  Rs  assuming 
a  straight  line  of  slope  Rs_1*2 
through  APS  =  25.  5  mbar  at  Rs 
=  9000  ft/kt*/3. 


Rs  (ft/kt1/3) 

AP  (mbar) 
s 

10 

89,  500 

30 

24,  000 

80 

7,  400 

100 

5,650 

300 

1,  510 

400 

1,  070 

500 

816 

600 

658 

700 

546 

1',  000 

357 

1,  500 

219 

2,000 

155 

3,  000 

95.4 

5,  000 

51.6 

6,  000 

41.5 

8,  000 

29.4 

9,  000 

25.5 

10.  000 

22.5 

15,  000 

13.8 

20.  000 

9.8 

50,  000 

3.26 

Figures  2,  3,  and  4  are  usedfor  this  pur¬ 
pose.  Enter  the  appropriate  figure  at  dob 
on  the  horizontal  axis,  go  up  to  the  appro¬ 
priate  curve  or  line,  and  read  fmax  from 
the  vertical  axis.  The  correct  figures  to 
use  are  summarized  below. 

Ground-Shock  f  _  Values  —  (deter- 

~  max  r /« 

mined  only  for  dob  >  60  ft/kt1'  ;  ground- 
shock -induced  airblast  is  not  important  for 
shallower  detonations) 

a.  Well-stemmed  nuclear  detonations 
in  dry  high-strength  rock,  60  <  dob  <  250 


to  300  ft/kt*/3,  0. 05  kt  <  W  <  5  kt.  Use 
the  nuclear  detonations  dry  rock  line 
(lower  line)  in  Fig.  2. 

b.  Nuclear  detonations,  nitromethane, 
TNT,  and  similar  chemical  explosive  det¬ 
onations,  in  alluvium  or  soil,  60  <  dob 

<  250  to  300  ft/kt1/3,  0.  001  kt  <  W  <  1.  0 
kt.  Use  the  alluvium  curve  (center  dashed 
curve)  in  Fig.  2. 

c.  Nitromethane,  TNT,  and  similar 
chemical  explosive  detonations,  in  dry 
high-strength  rock,  60  <  dob  <  250  to  300 
ft/kt1/3,  0.001  kt<  W<  1.0  kt.  Use  the 
chemical  explosives  in  basalt  or  rhyolite 
line  (upper  line)  in  Fig.  2. 

d.  Nitromethane  and  similar  chemical 
explosive  detonations  in  saturated  clay 
shale  and  saturated  weak  media,  at  opti¬ 
mum  depth  (dob  =  170  to  200 ft/kt1/3),  all 

yields.  Use  ground-shock  f  =  0.  06. 

max 

Gas-vent  airblast  is  negligible  for  these 
events.  The  airblast  from  detonations  in 
saturated  media  is  more  accurately  pre¬ 
dicted  by  the  theoretical  method  (discussed 
in  the  next  section). 

e.  Aluminized  ammonium  nitrate  slur¬ 
ry  and  ammonium  nitrate  fuel  oil  (ANFO) 
detonations  in  sandstone  or  weak  rock, 

160  <  dob  <  300  ft/kt1  ^3  0.001<  W<  0.1  kt. 
Use  the  ground-shock  f  line  in  Fig.  4. 
Airblast  from  detonations  of  this  type 
shallower  than  dob  =  160  ft/kt1/3  cannot  be 
accurately  predicted  due  to  lack  of  data, 
but  the  gas-vent  will  be  strongly  dominant, 
and  the  airblast  may  approach  that  ex¬ 
pected  from  a  surface  burst  of  the  same 
yield. 

Gas-Vent  f _ Values  — 

**  max 

a.  Nuclear  detonations  in  moist  or  wet 
alluvium,  moist  or  wet  rock,  end  all  moist 
or  weak  media,  15  <  dob  <  160  ft/kt1/3. 
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Maximum  transmission  factor 


Scaled  depth  of  burst  —  ft/kt^ 


Fig.  2.  Ground-shock  airblast  maximum  transmission  factor  fmax  fo1’  single-charge 
events  as  a  function  of  scaled  depth  of  burst  (for  TNT  and  similar  chemical 
explosives  in  basalt  and  rhyolite,  for  TNT  and  similar  chemical  explosives 
and  nuclear  explosives  in  alluvium,  and  for  nuclear  explosives  in  dry  high- 
strength  rock). 
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Maximum  transmission  factor 


Scaled  depth  of  bunt  —  ft  /let  ' 

Ground-shock  and  gas-vent  maximum  transmission  factors  fmax  for 
stemmed  single-charge  events  as  a  function  of  scaled  depth  of  burst  (for 
aluminized  ammonium  nitrate  slurry  explosive  in  weak  rock). 


0. 1  <  W  <  100  kt.  Use  the  solid  upper 
curve  in  Fig.  3. 

b.  Nitromethane,  TNT,  and  similar 

chemical  explosive  detonations  in  alluvium, 

soil,  and  high-  or  intermediate -strength 

rock,  15  <  dob  <  250to300ft/kt1/3,  0.001 

kt  <  W  <  1.0  kt.  Use  the  solid  upper  curve 

in  Fig.  3.  For  slightly  more  optimistic 

(lower)  estimates  of  f _ at  dob  >  160  ft/ 

.  I  o  max  ' 

ktx'  ,  use  the  dashed  line  in  Fig.  3. 

c.  Well-stemmed  nuclear  detonations 

in  dry'  high-strength  rock,  125  <  dob  < 
250  ft/kt1'3,  0.05  kt  <  W  <  5  kt.  Use  the 
straight  line  labeled  "f  (gas -vent)  for 

well-stemmed  nuclear  detonations  smaller 
than  5  kt  in  dry  high-strength  rock"  shown 
in  Fig.  3.  This  line  lies  well  below  the 
curve  (Cases  a  and  b,  above),  and  gas  - 
vent  airblast  is  usually  quite  small  for 
these  events.  Ground-shock-induced  air- 
blast  is  generally  dominant. 

d.  Aluminized  ammonium  nitrate  slur¬ 
ry  and  ANFO  detonations  in  sandstone  or 
weak  rock,  160  <  dob  <  30C  ft/kt1  , 

0.  001  <  W  <  0.  1  kt.  Use  either  of  the  two 
suggested  gas-vent  fmax  curves  (solid  or 
dashed)  in  Fig.  4.  The  upper  (solid)  curve 
will  give  a  somewhat  more  pessimistic 
(higher)  estimate  of  the  gas-vent  overpres  - 
sures.  Airblast  from  detonations  of  this 
type  shallower  than  dob  =  160  ft/kt1/3 
cannot  be  accurately  predicted  due  to 
lack  of  data,  but  gas-vent  overpressure 
will  be  strongly  dominant,  and  the  air- 
blast  may  approach  that  expected  from  a 
surface  burst  of  the  same  yield . 

The  fmajc  Values  for  Unstemmed 
Detonations — Buried  detonations  which 
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Moisture  content  less  than  or  equal  to 
1.0%  by  weight. 


are  unstemmed  produce  greater  airblast 

overpressures  than  do  normally  stemmed 

events  at  the  same  dob.  Unstemmed 

events  are  predicted  using  the  normal 

f _ method  and  the  f  values  given 

max  max 

below.  No  distinction  is  drawn  between 

ground-shock  and  gas-vent  airblast  for 

most  of  these  predictions,  and  the  fmax 

values  given  here  predict  the  greatest 

expected  airblast  overpressures.  (Note 

that  all  f  values  discussed  in  this 
max 

paragraph  are  based  on  experiments  with 
1-ton  aluminized  ammonium  nitrate 
slurry  chr.rges  in  weak  rock;  these  val¬ 
ues  may  not  produce  accurate  predictions 
for  widely  differing  yields,  explosive 
types,  or  media.)  To  predict  for  an  un¬ 
stemmed  chemical  explosive  detonation 
at  40  <  dob  <  360  ft/kt1/3,  with  a  cylindri¬ 
cal  open  shaft  to  the  surface  30  ft/kt1/3 

in  diameter,  use  f  »  0.99.  To  predict 
max 

for  an  unstemmed  chemical  explosive 

detonation  at  dob  >  200  ft/kt1/3,  with  a 

cylindrical  open  shaft  to  the  surface 

12  ft/kt1/3  in  diameter,  use  f  „  at  0,27. 

max 

To  predict  for  an  unstemmed  chemical 
explosive  detonation  at  dob  >  200  ft/kt1/3, 
with  an  open  shaft  to  the  surface  much 
less  than  12  ft/kt1/3  in  diameter:  because 
the  narrow  shaft  will  collapse  and  seal 
itself  at  early  time,  the  event  will  simu¬ 
late  a  stemmed  detonation;  therefore, 
predict  the  normal  ground-shock-induced 
and  gas -vent  airblast  for  a  stemmed 
detonation  at  the  same  dob. 

Step  4 

After  obtaining  the  ground- shock 

f  and  the  gas- vent  f _ values, 

max  max 

make  the  airblast  predictions  and 
plot  them  on  a  sheet  of  log- log  graph 
paper. 
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First,  predict  the  ground-shock-induced 

airblast  (detonations  deeper  than  dob 

=  60  ft/kt*^  only)  as  follows:  In  Step  2, 

several  overpressure  points  (R,  AP)were 

calculated  for  the  detonation.  Multiply 

each  overpressure  AP  by  the  ground-shock 

transmission  factor,  f  (ground- shock). 

max 

This  calculation  gives  the  predicted  ground - 
shock-induced  overpressure  points:  (R, 
AP  (ground-shock)). 

Next,  predict  the  gas-vent-induced 
airblast:  again,  use  the  overpressure 
points  (R,  AP)  calculated  in  Step  2.  Mul¬ 
tiply  each  overpressure  AP  by  the  gas- 
vent  transmission  factor,  f  (gas-vent). 

niaX 

This  calculation  gives  the  predicted  gas- 
vent- induced  overpressure  points:  (R, 

AP  (gas-vent)). 

The  overpressures  listed  in  Table  1 
are  in  millibars  (mbar).  Therefore,  the 
predicted  overpressures  will  also  be  in 
mbar.  If  desired,  convert  all  predicted 
overpressures  to  psi  by  dividing  them  by 
69.0: 

AP  (in  psi)  =  AP  (in  mbar)/69.  0. 

Step  5 

Plot  the  predicted  overpressures  on 
log- log  graph  paper  (predicted  overpres¬ 
sures  as  a  function  of  range,  R,  from 
SGZ).  First,  plot  the  ground -shoe k- 
induced  overpressures-  (R,  AP  (ground - 

shock)).  All  of  these  points  should  lie  on 

- 1  2 

a  straight  line  (of  slope  R  '  ).  Draw  a 

straight  line  through  them.  N^xt,  plot 

the  gas-vent-induced  overpressures:  (R, 

AP  (gas- vent  1).  Ail  of  these  points  should 

lie  on  another  straight  line  (also  of  slope 
- 1  2 

R  *  ).  Draw  a  straight  line  through 
them.  The  higher  of  the  two  straight  lines 
determines  the  peak  or  dominant  airblast 


overpressure  pulse  as  a  function  of  range. 
The  peak  overpressure  at  any  range  of 
interest  may  be  easily  predicted  by  means 
of  the  higher  line. 

Note  that  the  empirical  method  gives 

overpressures  at  ground-level  only.  AP 

is  predicted  as  a  function  of  R,  range 

along  the  ground  surface  from  SGZ.  Note 

also  that  the  predictions  are  not  valid  at 
1  2 

very  close  ranges.  '  The  overpressures 

may  be  somewhat  overpredicted  at  all 

ranges  such  that  the  original  Rg  (from 

Table  1,  before  scaling)  is  R  <  3000  ft/ 
1/3  s' 

kt  '  (gas-vent-induced  airblast)  or  R„ 

I/O  S 

<  600  ft/kt  '  (ground-shock-induced  air- 
blast).  In  other  words,  the  true  over¬ 
pressures  very  close  to  SGZ  will  be  less 
than  the  predicted  values.  A  more  thor¬ 
ough  review  of  close-range  airblast  will 
be  found  in  Ref.  1. 

PREDICTION  PROCEDURE  FOR  ROWS 
AND  ARRAYS  OF  CHARGES 

Airblast  overpressures  from  rows  or 
arrays  of  buried  charges  may  be  approx¬ 
imately  predicted  by  means  of  the  methods 
discussed  in  the  following  paragraphs.  The 
basic  procedure  is  to  predict  ground-shock- 
induced  airblast  and  gas-vent-induced  air- 
blast  from  an  "average"  single  charge  in 
the  row  or  array  (using  either  the  empiri¬ 
cal  or  the  theoretical  prediction  method). 
The  predicted  overpressures  are  then 
multiplied  by  appropriate  reinforcement 
factors.  A  few  important  differences 
between  singie-charge  and  multiple-charge 
predictions  should  be  mentioned.  First, 
all  multiple-charge  events  are  assumed 
to  be  detonated  simultaneously  unless 
otherwise  stated.  They  are  also  assumed 
to  be  emplaced  at  or  near  optimum 
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ajf  t 

intercharge  spacing  for  excavation  of  a 
continuous  crater  of  maximum  volume) 
unless  otherwise  stated.  The  range  from  a 
row  or  array  is  measured  horizontally  along 
the  ground  surface  as  usual,  but  SGZ  is 
considered  to  lie  at  the  geometric  center 
of  the  row  or  array.  Finally,  all  over¬ 
pressures  from  row  charges  are  predicted 
in  two  different  directions:  perpendicular 
to  the  row  axis  (1)  and  off  the  end  of  the 
row  (|| ).  The  reinforcement  factors  are 
always  different  for  these  two  directions, 
and  the  predicted  overpressure  at  a  given 
range  is  always  higher  perpendicular  to 
the  row  than  off  the  end.  The  overpressures 
at  directions  between  the  1  direction  and 
the  ||  direction  will  lie  between  the  pre¬ 
dicted  overpressures  for  the  two  directions. 

The  procedure  for  determining  the  re¬ 
inforcement  factors  and  predicting  air  - 
blast  from  a  row-charge  or  array  is 
summarized  below: 

Step  1 

Calculate  the  dob  of  each  charge  in  the 
row  or  array: 

dob  =  Unft/kt1/3) 

W1'3 

DOB  =  Depth  of  burst  to  center  of 
charge  (in  ft) 

W  =  Yield  of  charge  (in  kt) 

Take  the  average  of  the  scaled  depths  of 
all  the  charges.  This  average  is  consid¬ 
ered  to  be  the  mean  dob  for  the  entire 
row  or  array. 


Optimum  spacing  of  charges  is  approx¬ 
imately  equal  to  the  crater  radius  for  a 
single  charge  having  the  weight  of  one 
row-charge  member. 


Step  2 

Calculate  the  average  yield  of  the 
charges  in  the  row  or  array.  Sum  the 
yields  of  all  the  charges,  and  divide  the 
sum  by  the  total  number  of  charges,  n . 
This  average  is  considered  to  be  the  mean 
yield  per  charge  or  the  individual  charge 
yield  for  the  entire  row  or  array. 

Step  3 

Perform  a  normal  single -charge  air- 
blast  prediction  for  the  "average"  charge 
in  the  ro w  or  array.  Use  the  mean  dob 
calculated  above  for  dob,  and  the  mean 
yield  per  charge  for  W.  Predictions  may 
be  performed  by  the  Prediction  Procedure 
for  Single- Charge  Subsurface  Bursts,  or 
by  the  Theoretical  Prediction  Procedure, 
as  appropriate  for  the  particular  case  in 
question.  Begin  with  Step  2  of  either  pro¬ 
cedure. 

Step  4 

After  completing  the  single  -  charge 
predictions,  multiply  the  predicted  over¬ 
pressures  at  all  ranges  by  the  appropriate 
reinforcement  factors,  n  ,  where  n  is  the 
number  of  charges  in  the  row  or  array. 
Since  row- charge  reinforcement  factors 
for  the  1  direction  are  always  different 
from  the  factors  for  the  |]  direction,  sep¬ 
arate  predictions  must  be  made  for  these 
two  directions.  Also,  the  reinforcement 
factors  are  sometimes  different  for 
ground-shock-induced  overpressures  and 
gas-vent-induced  overpressures.  Be  sure 
to  use  the  proper  reinforcement  factor  in 
each  case.  The  following  is  a  summary 
of  the  reinforcement  factors. 

For  well-stemmed  nuclear  or  chemical 
explosive  row-charge  cratering  detonations 
of  large  individual  chargeyield  (W  >  0,  01  kt) 
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in  high-  or  intermediate- strength  rock: 


AP  (l  to  row)  =  n  '  AP  (single  charge  at 

same  range,  for 
ground-shock 
and  gas -vent 
overpressures) 

AP  (||  to  row)  =  n8*^5AP  (single  charge  at 

same  range,  for 
ground- shock 
overpressures) 

n  fi 

AP  (||  to  row)  =  n  *  AP  (single  charge  at 

same  range,  for 
gas-  vent  over¬ 
pressures). 

For  nuclear  or  chemical  explosive  row- 
charge  cratering  detonations  of  moderate 
to  large  individual  charge  yield  (W  >  0.  001 
kt)  at  or  below  optimum  depth  in  satura¬ 
ted  weak  media  such  as  clay  shale: 

0  7 

AP  (i  to  row)  =  n  *  AP  (single  charge  at 

same  range,  for 
ground- shock 
overpressures) 

AP  (||  to  row)  =  n8,41AP  (single  charge  at 

same  range,  for 
ground- shock 
overpressures). 

Ground-shock-induced  airblast  is  strongly 
dominant  for  both  row-  and  single-charge 
detonations  at  or  below  optimum  depth  in 
saturated  weak  media. 

For  chemical  explosive  row- charge 
cratering  detonations  of  small  individual 
charge  yield  (W  «  0.001  kt): 


0  9 

AP  (i  to  row)  =  n  *  AP  (single  charge  at 

same  range,  for 
both  ground- 
shock  and  gas- 
vent  over¬ 
pressures) 

AP  (||  to  row)  =  n8,  4  AP  (single  charge  at 

same  range,  for 
ground- shock 
overpressures) 


A  n 

AP  (||  to  row)  =  n  *  AP  (single  charge  at 

same  range,  for 
gas -vent  over- 
pressur'  ). 

For  all  square  array  detonations,  where 
n  is  the  total  number  of  charges  in  the  ar¬ 
ray  (n  =  4  or  5): 


AP  (for  array)  =  n0,65  to  n8,8  AP 

(single  charge  at  same  range,  for 
ground-shock  overpressures) 

AP  (for  array)  =  n8*8  to  n8*^  AP 

(single  charge  at  same  range,  for 
gas-vent  overpressures). 


The  above  relations  apply  in  all  directions 
from  the  geometric  center  of  the  array. 

For  chemical  explosive  row-  charge 
cratering  detonations  of  moderate  individ¬ 
ual  charge  yield  (W  =  0.  001  to  0.  01  kt)  in 
weak  or  moist  unsaturated  media  (such  as 
sandstone,  weak  rock,  or  soil),  inter  - 
charge  spacing  *=  160  to  350  ft/kt1/3: 


0  75 

AP  (l  to  row)  =  n  *  AP  (single  charge  at 

same  range,  for 
ground -shock 
overpressures) 

AP  ()|  to  row)  =  n°*  5  AP  (single  charge  at 

same  range,  for 
ground -shock 
overpressure  >) 

0  6 

AP  (l  to  row)  =  n  *  AP  (single  charge  at 

same  range,  for 
gas -vent  over¬ 
pressures) 


AP  (||  to  row)  =  n°*  8AP  (single  charge  at 

same  range,  for 
gas-vent  over¬ 
pressures). 


Recent  data  obtained  from  chemical  ex¬ 
plosive  detonations  of  moderate  individual 
charge  yield  in  weak  media  make  it 
possible  to  predict  reinforcement  factors 
ana  airblast  from  single-row  detonations 
with  delays  between  charges,  and  from 
simultaneous  and  delayed  double-row 


detonations.  The  prediction  procedures 
are  as  follows: 

Ground- Shock- Induced  Overpressures 
for  Inter  charge  Delays  (single  rows)  — 
Perpendicular  to  the  row,  intercharge  de¬ 
lays  greater  than  0.25  sec/kt1/**  (delays 
greater  than  0.025  sec  for  one-ton  charges): 

AP  (l  to  row)  =  aP  (single  charge  at 
same  range) — no 
reinforcement 
occurs. 

Off  the  starting  or  initiated  end  of  the  row 

(end  at  which  the  detonation  sequence  is 

1/3 

begun),  delays  greater  than  0.25  sec/kt  : 

AP  (||  to  row)  =  AP  (single  charge  at 
same  range). 

Off  the  final  or  concluding  end  of  the  row: 

Delay  =  0.25  S<T :  AP  (||  to  row 
kt1'*3 

0  5 

=  n  *  AP  (single  charge  at 
same  range). 

sec 

Delav  greater  than  0.  5  — 7-7-5-  : 

kt1/d 

AP  (||  to  row)  =  AP  (single 

charge  at  same  range). 

Gas-Vent-Induced  Overpressures  for 
Intercharge  Delays  (single  rows) —  Per¬ 
pendicular  to  the  row,  delays  greater  than 
0.  25  sec/kt1^: 

AP  (1  to  row)  =  ^  to  n®*  AP 

(sir  de  charge  at  same  range)  — 
estimate  only. 

Off  both  ends  of  the  row,  all  delays: 

AP  (||  to  row)  =  n0,  45  AP  (single 

charge  at  same  range). 

These  results  should  apply  most  reli¬ 
ably  to  events  near  optimum  depth  (dob 


*=  200  ft/kt1/  ^)  in  weak  media,  with  a. 
intercharge  spacing  =  250  ft/kt1^®. 

Double  Row  Events  —  For  optimum 
interrow  separation  (*350  ft/kt1^),  rein¬ 
forcement  factors  and  airblast  are  ap¬ 
proximately  the  same  as  for  a  single  row 
containing  the  same  number  of  charges 
as  both  of  the  double  rows  (n  =  total 
number  of  charges  in  both  rows).  The 
SGZ  is  considered  to  lie  at  the  geometric 
center  of  both  rows,  and  predictions  are 
made  as  discussed  above.  At  wide  inter¬ 
row  separations  (greater  than  1100  ft/ 
kt1/**)  or  long-interrow  delay  times  for 
optimum  separation  (delay  time  between 
rows  greater  than  1.  2  sec/kt1/**),  the  air- 
blast  pulses  from  the  two  rows  will  not 
normally  combine  or  reinforce.  Thus  , 
airblast  may  be  predicted  for  each  row 
separately,  as  though  that  row  were  being 
detonated  alone.  The  respective  SGZ's 
of  the  two  rows  are  considered  to  lie  at 
the  geometric  center  of  each  row.  Pre¬ 
dictions  are  made  as  discussed  above,  for 
each  of  the  two  rows  separately.  This  pre¬ 
diction  technique  may  prove  incorrect  for 
a  direction  in  which  the  distances  or  delay 
times  between  any  of  the  various  pulses 
from  the  two  rows  is  insufficient  to  pre¬ 
vent  overlap  of  the  pulses,  or  for  very 
long  rows  with  spatially  extended  over¬ 
pressure  pulses  (total  row  length  ^un¬ 
sealed  delay  time  X  local  sonic  velocity). 
There  is  currently  no  good  method  for 


For  example,  pulses  tend  to  overlap 
along  the  exact  center  line  direction  off 
the  end  of  a  simultaneously  detonated  dou¬ 
ble  row;  overlap  may  also  occur  perpen¬ 
dicular  to  a  delayed  double  row,  in  the 
perpendicular  direction  closer  to  the 
later-detonated  of  the  two  rows. 
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predicting  airblast  from  yery  long  double 
rows  consisting  of  many  charges.  \ 

Airblast  irom  Row  Charges  on  Sloping 
Ground-— For  row- charge  events  on  mod¬ 
erate  hillside  slopes,  there  is  no  signifi- 
cant'difference  between  peak  overpressures 
in  the  uphill  and  downhill  directions  from 
the  row  (the  uphill  results  appear  to  be 
very  slightly  higher  in  some  cases).  Pre¬ 
dictions  are  performed  as  usual  (for  an 
event  on  level  ground),  except  that  ranges 
are  measured  along  the  sloping  ground 
surface  from  the  geometric  center  of  the 


After  multiplying  the  overpressures 
t^y  the  appropriate  reinforcement  factors, 
plot  the  reinforced  overpressures  as  a 
function  of  range  on  log-log  graph  paper. 
Plot  both  the  ground- shock- induced  and 
gas -vent- induced  overpressures.  Each) 
set  of  overpressures  will  define  a  straight1 
line.  The  higher  of  the  two  lines  at  a 
given  ang^  will  determine  the  ^iominant 
airblast  pulse  and  peak  overpressures . 

In  the  case  of  row-charge  events,  sepa¬ 
rate  airblast  predictions  are  made  for  the 
1  and  ||  directions.  ■  Both  the  ground- 
shock- induced  and  the  gas -vent- induced 
overpressures  should  be  plotted  for  each 
direction.  Predictions  for  the  two  direc¬ 
tions  should  be  plotted  separately.  The 
higher  of  the  two  lines  at  a  given  range 
for  a  given  direction  will  determine  the 

.  I 

dominant  airblast  pulse  and  the  peak  over¬ 
pressures  in  that  direction. 

PREDICTION  PROCEDURE  FOR  SURFACE 
BURSTS  AND  FREE- AIR  BURSTS 

The  ^irblast  overpressures  from!  near- 
surface ’bursts,  buried  detonations  at 


scaled  depths  shallower  than  15  ft/kt1^, 
and  free-air  bursty  are  also  predicted  by', 

;  a  scaling  procedure.  However)  the  stand¬ 
ard  l^ine  point?  listed  in  Table  1  cannot 
be  used  because  a  straight  line  fit  is  in¬ 
adequate  for  these  events.  Instead  Figs. 

5  and  6  are  used.  These  figures  are  based 
on  the  best  available  data  from  recent  ex- 
perimenta'J  chemical  explosive  tests,  and 
on  theoretical  shock  calculations.  * 

To  predict  airblast  from  near-surface 
bursts,  buried  detonations  at  dob  <15 
ft/kt1^,  or  free-air  bursts,  read  the 
overpressures,  APg,  at  several  ranges, 

Rg  J  from  the  appropriate  curve  or  set  of 
points  (in  Fig.  5  or  Fig.  6).  Scale  these 
overpressures  and  ranges  to  the  yield, 

W,  and  ainbient  pressure,  P.;m5ient,  of 
the  detonation  to  be  predicted,  using  Eqs.i 
(2)  and  (3)  —  see  the  section:  Prediction 
Procedure  for  Single- Charge  Subsurface 
Bursts.  The  calculated  points  (R,  AP) 
give  the  predicted  airblast  for  the  detona¬ 
tion.  These  points  may  be  plotted  on  log- 
log  graph  paper,  and  a  smooth  curve  may) 
be  drawn  through  them.  The  points  will 
not  lie  on  a  straight  line  in  th^s  case . 
Several  points  at  various  ranged  should 
be  calculated  and  plotted  in  order  to  de¬ 
fine  the  curve  properly. 

The  following  is  a  discussion  of  the 
correct  curves  to  use  for  predicting  vari¬ 
ous  cases.  For  TNT  and  similar  chemi¬ 
cal  explosive  detonations  near  the  surface 
and  at  dob's  shallower  than  15  ft/kt1^, 
0.1001  <  W  <  1.0  kt,  use  the  upper  curve 
in  Fig.  5  or  the  set  of  points  just  below 
the  upper  curve  in'  Fig.  5.  The  points 
provide  a  slightly  lower  (lesfe  pessimistic) 
estimate  of  the  predicted  overpressures. 

For  ANFO  detonation;  near  the  surface 
and  at  dob's  shallower  than  15  ft/kt5^  , 


.  V 
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\ 
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mbar 


Fig.  5.  Overpressure  as  a  function  of  range  for  nuclear  and  TNT  chemical  explosive 
surface-burst  events  (yield  =  1.0  kt,  ambient  pressure  =  1000  mbar,  and  for 
all  free-air  bqrst  events  (yield  =  1.0  kt,  ambient  pressure  =  1000  mbar). 
Circles  represent  minimum  observed  data  points  for  multi-ton  TNT  surface 
'  bursts. 

\ 

1 


0.  001  <  W  <  1.0  kt,  use  the  plotted  points  . 
(squares  and  stars  only)  in  Fig.  6.  Beyond 
the  range  bf  the  plotted  points,  use  the 
dashed  curve  in  Fig.  6.  t 

For  aluminized  ammoriium  nitrate 
slurry  detonations  near  the  surface  and 
at  dob's  shallower  than  15  ft/kt * 0.001 
<  W<  1.0  kt,  use  the  solid  and  dashed 
curve  in  Fig.  6. 

For  free-air-burst  detonations  (events 
in  a  homogeneous  atmosphere  with  nb 


nearby  ground  surface)  of  all  types  and 
yields,  use  the  lower  solid  curve^in  Fig.  5, 
If  the  event  is  chemical  explosive  multi¬ 
ply  the  true  yield,  W  (in  kt),  by  a  factor 
of  2.0  before  substituting  it  into  Eq.  (3); 
if  the  event  is  nuclear  explosive,  use 
the  true  yield,  W.  The  free-air-burst 
curve  shown  in  Fig.  5  is  based  on  a 
theoretical  calculation  for  nuclear 
free- air  bursts,  known  as  the  IBM 
Problem  M. 


Fig.  6.  Overpressure  as  a  function  of  range  for  typical  ammonium  nitrate  fuel  oil  and 
aluminized  ammonium  nitrate  slurry  surface-burst  events  (yield  =  1.0  kt, 
ambient  pressure  =  1000  mbar).  TNT  surface-burst  event  curve  is  included 
for  comparison. 


SAMPLE  PROBLEM  FOR  SINGLE- 
CHARGE  SUBSURFACE  BURST 


Go  to  Table  1  and  select  a  few  appro¬ 
priate  points  at  ranges  of  interest: 


Predict  the  ground- shock- induced  and 
gas-vent-induced  airbiast  overpressures 
for  a  20-ton  (0.02  kt)  TNT  detonation  at 
a  DOB  of  34  ft  in  alluvium.  The  ambient 
pressure  at  the  project  location  is  Pamt,jent 
=  868  mbar.  First,  find  the  scaled  depth 
of  burst,  dob: 


dob 


DOB  .  34  ft 

W1^3  (Q.02)1^3 


=  125ft,ikt1//3. 


R  (ft)  APq  (mbar) 

500  816 

1000  357 

3000  95.4 

Scale  these  points  back  to  a  20- ton  ex 
periment  at  Pambien,  ■  868  mbar: 

&  =  ^a  (Ml)  '  °- 868 
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R  --  R 


s 


/O,  02 
\  1.0 


iooo\ 1//3 

868  ) 


--  0.  285  Rd. 

s 


Tabulate  the  values  of  R  and  aP: 


R  (ft) 

AP  (mbar) 

142 

709  =  10. 28  psi 

285 

310  =4. 49  psi 

855 

82.8  =  1.  20  psi 

The  experiment  is  in  alluvium.  Use 
the  dashed  (alluvium)  curve  in  Fig.  2. 
Entering  the  figure  at  dob  =  125  ft/kt1/3, 
find  f  (ground-shock)  -  0.  03.  Using 
the  upper  curve  in  Fig.  3  and  entering  the 
figure  at  dob  *  125  ft/kt1/3,  find  f  „ 

max 

(gas-vent)  =  0.  17. 

For  ground-shock  overpressure  pre  - 
dictions,  multiply  the  overpressures 
(either  in  psi  or  mbar,  as  desired)  by 


* 


To  convert  to  psi,  divide  by  69. 


f  (ground- shock)  =  0.03  and  tabulate: 
max 


R  (ft) 

AP  ground-shock  (psi) 

142 

0.308 

285 

0.  135 

855 

0.036 

For  gas- 

vent  overpressure  predictions, 

multiply  the  overpressures  (either  in  psi 

or  mbar,  as  desired)  by  f  (gas-vent) 

J  max 

=  0. 17  and  tabulate: 

R  (ft) 

AP  gas-vent  (psi) 

142 

1.75 

285 

0.764 

855 

0.204 

These  gas -vent  and  ground -shock  pre¬ 
dictions  are  plotted  in  Fig.  7.  They  defin» 
two  straight  lines,  which  are  the  predicted 
gas-vent  and  ground-shock  overpressures 
Gas-vent- induced  airblast  overpressures 
are  dominant  in  this  case. 


Theoretical  Prediction  Method  (for  Ground-Shock-lnduced  Airblast) 


THEORETICAL  PREDICTION  PROCEDURE 

The  theoretical  method  predicts  ground- 
shock  -  induced  airblast  overpressures 
from  most  types  of  buried  single-charge 
detonations  at  60  <  dob  <  700  ft/kt1/3.  It 
does  not  predict  for  row -charge  events, 
but  it  can  be  used  to  perform  the  ground - 
shock-induced  single-charge  predictions 
which  are  then  multiplied  by  the  ground- 
shock  reinforcement  factors  to  predict 
ground-shock-induced  airblast  from  row- 
charge  events  (see  the  previous  section 
under  Prediction  Procedure  for  Rows  and 
Arrays  of  Charges). 

The  theoretical  method  does  not  utilize 
scaling  relations.  Instead,  a  peak  local 
overpressure,  APq,  directly  above  SGZ 


is  calculated  from  the  peak  vertical  spalls 
velocity  at  SGZ,  VQ.  This  local  overpres¬ 
sure  is  used  in  predicting  overpressures 
at  all  other  ranges.  APq  is  determined 
from  the  equation: 


^0  '  1013  mbar  (p0c0)MV0‘ 


(4) 


where 

AP_  =  ground -shock- induced 
u  peak  local  overpressure 
directly  above  SGZ  (in 
psi) 

P  hi  t  =  ambient  atmospheric  pres¬ 
sure  at  location  of  exper¬ 
iment  (in  mbar) 

p-c.  =  acoustic  impedance  of  air 
at  standard  sea*  level 
conditions 

--  0.01841  P8lY~t3ec 


-19- 


cQ  =  sonic  velocity  in  air  at 
u  sea  level 

=*  1087  ft/sec 

VQ  =  peak  vertical  spall  veloc¬ 
ity  of  ground  surface  at 
SuZ  (in  ft/ sec) 

M  =  Mach  number  of  air  shock 
at  SGZ  (M  *»  1. 0  for  low 
velocities,  V^;  see  below ). 

The  quantities  appearing  in  Eq.(4)  are 
discussed  in  detail  in  the  following  para¬ 
graphs.  If  APq  is  known,  the  surface 
level  overpressure,  AP,  at  any  range,  R, 
from  SGZ  can  be  calculated  from  the 
equation: 

^  '  “o  “(tt)'  (5) 

where 

AP  =  peak  ground-shock-induced 
overpressure  (in  psi) 

DOB  =  depth  of  burst  to  center  of 
charge  (in  ft) 

R  =  range  along  ground  surface 
from  SGZ 

a  =  a  constant  for  all  detonations 
in  a  given  medium  (dimension¬ 
less). 

The  theoretical  method  is  also  capable 
of  predicting  ground- shock- induced  air  - 
blast  at  elevated  locations  (above  ground 
surface  level).  Assume  that  an  elevated 
point  is  located  at  height  h  above  the 
ground  surface  plane,  and  at  range  R 
from  SGZ  (this  is  a  true  range  or  slant 
range,  not  the  range  along  the  ground 
surface).  The  overpressure  AP  (elevated) 
in  psi  at  this  point  is  given  by: 

AP  (elevated)  =  APq  or 

x[cos  75^-gy]’  (6) 


where:  R  =  slant  range  from  SGZ  to  the 
elevated  point  (in  ft),  and  6  is  determined 
from: 

cos  6  =  .  (7) 

Equation (6)  is  not  valid  if  (90°  -  0)  >  7  5°  (i.  e. , 
for  elevated  points  located  almost  directly 
above  SGZ). 

In  order  to  apply  Eqs.  (5)  and  (6),  it 
is  necessary  to  know  the  value  of  the  con¬ 
stant  a.  The  value  of  a  is  controlled  by 
the  shape  of  the  vertical  velocity  field 
above  a  detonation  (in  other  words,  by  the 
vertical  velocity  profile  at  locations  away 
from  SGZ).  The  shape  of  the  vertical  ve¬ 
locity  field  is  the  same  for  all  single  - 

charge  experiments  at  all  dob*  s  in  a  given 
2 

medium.  The  shape  does  vary  some  - 

what  for  different  media.  More  precisely, 

the  peak  vertical  spall  velocities  above 

single- charge  detonations  at  locations 

other  than  SGZ  are  accurately  fitted  by  an 
2 

equation  of  the  form  : 


V  =  peak  vertical  spall  velocity  of 
ground  surface  at  slant  range 
S  from  shot  point 

Vn  *  peak  vertical  spall  velocity  at 
u  SGZ 

DOB  =  depth  of  burst  of  charge  (in  ft) 

S  =  slant  range  from  the  shot  point 
to  the  ground  surface  location 
(in  ft) 

m  =  a  velocity  profile  constant  for 
all  single-charge  detonations 
in  a  given  medium 

It  can  be  shown  that,  for  single* charge 
velocity  fields  of  this  type,  the  constant 
a  in  Eqs.  (5)  and  (6)  is  a  function  only  of 
m.  If  the  velocity  field  constant  m  for  a 
particular  medium  is  known,  the  value  of 
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a  can  be  determined  by  a  straightforward 
analytic  solution.  The  details  of  this  anal¬ 
ysis  are  discussed  in  Ref.  2.  The  analy¬ 
sis  has  been  carried  out  for  typical  values 
of  m,  and  the  resultant  a*s  are  listed  in 
Table  2.  The  first  column  of  the  table 
lists  the  value  of  m.  The  second  column 
gives  the  corresponding  value  of  a.  The 
final  column  lists  the  medium  to  which 
this  value  of  a  applies.  The  a1  s  listed  in 
Table  2  will  be  used  in  making  theoretical 
airblast  predictions. 

Given  the  value  of  a  for  a  particular 
medium,  it  is  possible  to  predict  ground- 
shock-induced  airblast.  The  required  in¬ 
put  data  for  the  prediction  include: 

(1)  Medium  and  approximate  moisture 
content 

(2)  Yield  W  of  the  explosive  charge 
(in  kt) 

(3)  Depth  of  burst  to  center  of  charge 
"DOB"  (in  ft) 

(4)  The  ambient  atmospheric  pressure 
at  the  location  of  the  project, 
Pambient  (in  mbar).  The  ambient 
pressure  may  be  estimated  from 
meteorological  data  or  from  the 
altitude  of  the  project  location  above 
sea  level:  (Pambient  58  1013  mbar 
at  sea  level,  Pambient  *  844  mbar 
at  5000  ft  above  sea  level). 


(5)  The  value  of  a  (from  Table  2). 

The  steps  in  the  prediction  procedure 
may  be  summarized  as  follows: 

Step  1.  Calculate  the  scaled  depth  of 
burst  dob  for  the  detonation: 

dob  =  (dob  in  ft/kt1/3) 

W1'3 

Step  2.  Determine  c  for  the  medium 
of  the  project  (Table  2). 

Step  3.  Determine  the  peak  vertical 
spall  velocity  at  SGZ,  VQ  (in  ft/sec). 
Peak  vertical  spall  velocities  are  plotted 
as  a  function  of  dob  in  Figs.  8  through  12. 
Enter  the  appropriate  figure  at  dob  on  the 
vertical  axis,  go  across  to  the  appropriate 
line,  and  read  Vq  (in  ft/sec)  from  the  hor¬ 
izontal  axis.  The  following  is  a  summary 
of  the  correct  figures  to  use: 

For  all  detonations  in  alluvium  and  soil, 
60  <  dob  <  300  ft/kt1/3,  use  Fig.  8. 

For  nuclear  detonations  in  dry  rhyolite, 
60  <  dob  <  300  ft/kt1/3,  use  Fig.  8. 

For  chemical  explosive  detonations  in 
dry  high-strength  rock  (basalt,  rhyolite, 
etc.),  60  <  dob  <  300  ft/kt1^3,  use  the 
upper  line  in  Fig.  9. 


Table  2.  Values  of  "a"  for  several  media  of  interest. 


Velocity  profile 
constant,  m 

a 

Medium 

2.0 

0.  50 

Maximum  value  for  contained  detonations  in 
any  medium 

2.  13 

0.46 

Underwater  detonations  and  underwater  det¬ 
onations  in  saturated  weak  media,  such  as 
coral 

2.  5 

0.3815 

Saturated  clay  shale  and  other  saturated 
weak  media 

3 

0.318 

Alluvium,  soil,  and  rhyolite  rock 

4 

0.250 

5 

0.212 

(Basalt  rock  ?) 

6 

0.  188 

Basalt  rock 

21- 


Scaled  depth  of  burst  —  ft/kt 


Fig.  8.  Peak  vertical  spall  velicity  as  a  function  of  scaled  depth  of  burst  for  all  events 
in  alluvium  and  for  nuclear  events  in  unsaturated  rhyolite. 
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Peak  vertical  spall  velocity  at  SGZ,  Vq  —  ft /sec 


Fig.  9.  Peak  vertical  spall  velocity  as  a  function  of  scaled  depth  of  burtt  for  chemical 
explosive  events  in  basalt  and  rhyolite,  and  for  nuclear  events  in  basalt. 
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Pig.  10.  Peak  vertical  spall  velocity  as  a  function  of  scaled  depth  of  burst  for  all 
events  in  weak  saturated  rock,  such  as  saturated  clay  shale. 


Peak  vertical  spall  velocity  at  SGZ,  Vq  —  ft/sec 

Fig.  11.  Peak  vertical  spall  velocity  as  a  function  of  scaled  depth  of  burst  for  events 
in  seawater  and  for  events  in  a  very  weak  saturated  medium  such  as  weak 
coral  overlain  by  seawater. 
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Fig.  12.  Peak  vertical  spall  velocity  as  a  function  of  scaled  depth  of  burst  for  all  coi 
tained  and  mounding  events  (scaled  depths  of  burst  greater  than  300  ft/kt */ 


For  nuclear  detonations  in  dry  basalt, 
60  <  dob  <  300  ft/kt1/3,  use  the  lower 
line  in  Fig.  9. 

For  all  detonations  in  weak  saturated 
media  (such  as  saturated  clay  shale),  60 

<  dob  <  300  ft/kt1/3,  use  Fig.  10.  Ground- 
shock-induced  airblast  is  strongly  domi¬ 
nant  for  all  such  events  deeper  than  dob 

=  170  ft/kt1/3. 

For  all  detonations  in  water,  80  <  dob 

<  700  ft/kt1/3,  use  the  upper  line  in  Fig. 
11.  Ground- shock- induced  airblast  (more 
correctly,  water-shock- induced  airblast) 
is  strongly  dominant  for  all  such  events 
deeper  than  dob  =  170  ft/kt1/3. 

For  all  detonations  in  very  weak  satu¬ 
rated  media  overlain  by  water  (such  as 
very  weak  saturated  submerged  coral),  80 

<  dob  <70  ft/kt1/3,  use  the  lower  line  in 

Fig.  11.  Ground- shock- induced  airblast 

(more  correctly,  water-shock-induced 

airblast)  is  strongly  dominant  for  all  such 

1  /  3 

events  deeper  than  dob  =  170  ft/kt  '  . 

For  all  contained  and  mounding  deto¬ 
nations  (normally,  events  at  700  >  dob 
§  300  ft/kt1/3),  use  Fig.  12  (lower  line  for 
'.•vents  in  alluvium  and  soil;  upper  line  for 
all  other  media).  Ground-shock-induced 
airblast  is  dominant  for  all  contained  or 
mounding  detonations. 

Step  4.  Calculate  (Vq/c),  where  c 
=  local  sonic  velocity  in  air.  (c  ~  1087 
ft/sec  for  standard  sea- level  condi¬ 
tions,  and  deviates  from  this  value  by 
only  a  few  percent  for  most  cases  of 
interest. ) 

Step  5.  Determine  the  Mach  number, 
M,  of  the  shock  (Ref.  2)  using  Fig.  13a 
(for  small  values  of  Vq/c)  or  Fig.  13b 
(large  values  of  Vq/c).  Enter  the  figure 
at  the  calculated  v  ilue  of  (Vq/c)  on  the 
vertical  axis,  go  across  to  the  curve,  and 


read  the  Mach  number,  M;  from  the  hori¬ 
zontal  axis. 

Step  6.  Using  the  above-determined 
values  of  VQ,  and  M  in  Eq.  4, 

calculate  the  SGZ  local  overpressure, 

APQ  (in  psi). 

Note:  PqCq  =  0.01841  PSA~?ec 

Step  7.  Substitute  the  predicted  value 
of  APq  into  Eq.  (5),  and  calculate  the 
overpressures,  AP  (in  psi),  at  several 
ranges  of  intt-est,  R. 

Step  8.  Plot  the  predicted  overpres¬ 
sures,  AP,  as  a  function  of  R  on  log-log 
graph  paper.  The  plotted  points  will  de¬ 
fine  a  straight  line  of  slope  R-1**1.  Draw 
a  line  through  the  predicted  points.  This 
line  may  be  used  to  predict  the  ground- 
shock-induced  overpressures  as  a  func¬ 
tion  of  range.  The  predictions  are  not 
valid  for  (R/'DOB)  <  1.0  (most  cratering 
detonations  in  unsaturated  media)  or  for 
(R/  DOB)  <  2.  0  (detonations  in  saturated 
media,  detonations  under  water  and  in 
media  overlain  by  water,  and  contained 
or  mounding  detonations).  At  these  ver> 
close  ranges,  APq  provides  an  approxi¬ 
mate  estimate  of  the  peak  ground-shock- 
induced  overpressure  near  the  ground 
surface. 

Step  9.  If  desired,  predict  the  ground- 
shock-induced  overpressures  AP  (elevated) 
at  any  elevated  locations  of  interest,  using 
Eqs.  (6)  and  (7).  Elevated  overpressure 
predictions  are  not  valid  if  (90°  -  6)  >  75°. 

SAMPLE  PROBLEM 

Predict  the  ground-shock-induced  air- 
blast  overpressures  from  a  20-ton  (0.02- 
kt)  single- charge  TNT  detonation  at  DOB 
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=  34  ft  in  a;  luvium.  The  ambient  atmos¬ 
pheric  pressure  at  the  experiment  loca- 

tionisPambienf=  867  mbar*  ' 

Yield  W  =  20-ton  chemical  explosive 
»  0. 02  kt 

Depth  of  burst,  DOB  =  34  ft 

P  ambient  =  867  mbar  (known  from 
meteorological  data  or  altitude) 

Scaled  depth  of  burst,  dob 

_  DOB  (iri  ft)  .  34 

[W  (in  kt)]  (0.02)T73 

«  =  125  ft/kt1/3 

o  =  0.318  (for  alluvium.  Table  2) 

Referring  to  Fig.  8  ^for  alluvium)  at 

dob  *  125  ft/kt*/3,  the1  predicted  SGZ 

peak  vertical  spall  velocity,  Vq,  is  138 

ft/sec.  The  local  sonic  velocity  c  =*  1087 

ft/s^c.  ‘  Therefore,  ^0  =  o.  127.  Entering 
1  c 

Fig.  13b  with  this  value,  the  Mach  num¬ 
ber  M  is  found  to  be  M  -  1.  079.  Now  it 
is-  possible  to  calculate  the  local  SGZ 
overpressure,  ,AP0,  from  Eq.  (4): 


”  V 


APo  *  refflifr  'Vo1  MVo  ■  2-347  ps1- 


Substituting  this  value  into  Eq.  (5), 

aP=APo0(OOB)  I 

=  (2.347)  (0.318)  (34)/R  =  25.4/R. 

^Sample  predicted  overpressures  from 
this  equation  are  tabulated  below: 

Distance  fropi  SGZ, 

_ R  (in  ft) _ AP  (psi) 


25 
50 
100 
300 
1000  \ 


1.02 
0.  51 
0.  25 
0.  085 
0.025 


,  l 

Figure  14  shows  the  predicted  line  drawn 

through  these  sample  points.  Its  slope  is 


Fig.  14.  Sample  problem:  ground-shock- 
induced  airblast  prediction, 
using  theoretical  method. 


\ 


\ 


R  *.  Do  not  extepd  predictions  inside 
R/DOB  =  1.0.  At  ranges  (R/DOB)  less^ 
than  1.0,  APq  provides  an  estimate  of  the. 
ground-shock-induced  peak  overpressure 
(APq  =  2.347  psi). 

The  theoretical  method  also  predicts 
elevated  overpressures.  An  imaginary 
problem  has  been  constructed  to  demon¬ 
strate  this  application.  Predict  the 
ground-shodk-indueed  peak  overpressure 
at  an  elevated  gage  100  ft  high,  at  true  1 
range  R  =  300  ft  from  SGZ: 

\ 

h  100  ft. 

Therefore,  the  angle  0  from  the  vertical 
to  the  gage  direction  (see  sketch)  is 
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AP  (elevated)  =  (2.347)  (0.318)(^) (1.06) 
AP  (elevated)  =  0.09  psi  at  elevated  gage. 

This  concludes  the  theoretical  predic¬ 
tion  of  ground -shock-induced  airblast. 
Gas -vent  airblast  must  still  be  predicted 
by  the  empirical  method,  as  discussed  in 
the  previous  section. 


Summary 

An  empirical  method  is  presented  for  of  the  two  methods  are  reviewed,  and 
predicting  both  ground-shock-induced  and  suggestions  for  applying  the  methods  in 
gas -vent-induced  airblast  overpressures  practice  are  set  forth.  Sample  problems 
from  buried  single-  and  row-charge  era-  are  included  to  clarify  the  prediction  pro- 

tering  detonations  in  a  variety  of  media.  cedures.  This  report  summarizes  state- 

A  semitheoretical  technique  for  predict-  of-the-art  airblast  prediction  techniques 
ing  ground-shock-induced  airblast  from  for  typical  underground  detonations,  but 
single-charge  events  is  also  discussed.  does  not  provide  detailed  supporting  data 

The  range  of  applicability  and  limitations  or  discussions. 
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